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IBM SYSTEM/360 
ASSEMBLY LANGUAGE 
INTERVAL ARITHMETIC SOFTWARE 


INTRODUCTION 

This document describes some computer software designed to perform interval 
arithmetic. An interval is defined as the set of all real numbers between two 
given numbers including or excluding one or both endpoints. Interval arithmetic 
consists of the various elementary arithmetic operations defined on the set of all 
intervals, such as interval addition, subtraction, union, etc. The theoretical 
aspects of interval arithmetic are described in reference (2). 

One of the main applications of interval arithmetic is in the area of error 
analysis of computer calculations. For example, it has been used successfully 
to compute bounds on sounding errors in the solution of linear algebraic 
systems^ 2) , error bounds in numerical solutions of ordinary differential 
equations/ 1)(3) as well as integral equations and boundary value problems. 

The described software should enable users to implement algorithms of the 
type described in these references efficiently on the IBM 360 system. 

The subroutines in this package were supplied to Goddard Space Flight 
Center by Dr. Ronald Tost, Institute fur Numerische Datenverarbeitung under 
an~exchange agreement with Dr. C. E. Velez, GSFC Code 552. Although the 
subroutines were coded in IBM 360 assembly language, they were designed so 
that calls may be made from FORTRAN programs. The sample problem in 
Section Vn and all of the examples in Sections I to VI are given in FORTRAN 
to illustrate how this package may be used in FORTRAN programs. The sub- 
routines were tested under varying conditions by Computing and Software, Inc., 
personnel . 

The following pages explain in detail the arithmetic operations performed 
by this package . 
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INTERVAL ARITHMETIC PACKAGE 


I. Language 

Assembler Language Code for IBM/360. 


II. Purpose 

This package is designed to perform arithmetic operations using floating 
point (64 bits) interval arithmetic. 


III. Method 

The interval members are defined as COMPLEX * 16, by the user. The basic 
arithmetic operations are performed along with some basic set theory applica- 
tions. The arithmetic functions are: addition, subtraction, multiplication, divi- 
sion, exponentiation. The absolute value of an interval and a negative interval 
can be created. Real and integer numbers can be converted into interval numbers. 
The intersection and union between two sets of interval numbers can be found. 

An illegal interval can be tested for. There are tests to determine if zero is 
contained in (e ) an interval or if one interval is contained in another. The sum of 
two sets of interval numbers can be found by means of a scalar product. The 
above operations are accomplished by the following functions: 

DLI, DRE, $IT$, $DT$, $ADD, $SUB, $MUL, $DIV, $POT, $MIN, $ABS, $SCH, $VER, 
IZER, ICOR, $ENT, $SUM, $INP. 


IV. Calling Sequence 


The interval numbers and functions are defined as COMPLEX * 16. 

In the following set: A is an INTERVAL (AL, AR) 

B is an INTERVAL (BL, BR) 

C is an INTERVAL (CL, CR) 

D: REAL * 8 
N: INTEGER * 4 


Functions: 

1. D = DLI(A) D: = AL Left Corner of A 
Ex. A = (2.5, 4.0) 

D = DLI(A) = 2.5 
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Right Corner of A 


2. D = 
Ex. 


3. A = 
Ex. 


4. A = 
Ex. 


5. A = 
Ex. 


6. A = 
Ex. 


7. A = 
Ex. 


8. A = 
Ex. 


DRE(A) D: = AR 
A = (2.5, 4.0) 

D = DRE(A) = 4.0 

$IT$(N) A: = (N, N) Produce an interval from an integer 
N = 5 

A = $IT$(N) = (5.0, 5.0) 

$DT$(D) A: = (D, D) Produce an interval from a real number 
D = 2.52 

A - $DT$(D) = (2.52, 2.52) 

A = $DT$ (.5DO) = (0.5D0, 0.5D0) 

$ADD (B, C) A: = B + C Addition 
B - (3.0, 4.0) 

C = (2.0, 2.0) 

A = $ADD (B, C) = (3.0 + 2.0, 4.0+ 2.0) 

= (5.0, 6.0) 

B = (-1.5, -0.5) 

C - (2.5, 4.5) 

A = $ADD (B, C) = (-1.5 + 2.5, -0.5 + 4.5) 

= (1.0, 4.0) 

$SUB(B,C) A: = B - C Subtraction 
B = (3.0, 4.0) 

C = (2.0, 2.0) 

A = $SUB (B, C) = (3.0 - 2.0, 4.0 - 2.0) 

= (1.0, 2.0) 

B = (3.0, 4.0) 

C = (-1.5, -2.5) 

A = $SUB (B, C) = (3.0 - (-1.5), 4.0-(-2.5)) 

= (4.5, 6.5) 


$MUL (B, C) A: = B * C Multiplication 
B = (3.0, 4.0) 

C = (2.0, 2.5) 

A = $MUL (B, C) = (MIN (6.0, 7.5, 8.0, 10.0), MAX (6.0, 7.5, 8.0, 10.0)) 

= (6.0, 10.0) 

$DIV (B, C) A: = B/C Division 
B = (3.0, 4.0) if 0 i (B, C) 

C = (1.0, 2.0) 

A = $DIV (B, C) = (3.0, 4.0) * (0.5, 1.0) 

= (MIN (1.5, 3.0, 2.0, 4.0), MAX (1.5, 3.0, 2.0, 4.0)) 

= (1.5, 4.0) 
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9. A = $DOT (B, N) A: = B**N Exponentiation 
Ex. B = (3.0, 4.0) 

N = 2 

A = $DOT (B, N) = (9.0, 16.0) 


10. A = $MIN (B) A: = - B Negative Interval 
B = (3.0, 4.0) 

A = $MIN(B) = (-4.0, -3.0) 


11. A = $ABS (B) 


A 

A 

A 


= B if BL > 0 Absolute Value 
= - B if BR <0 else 
= (0, MAX(BL, BR)) 


Examples: . 

(1) B = (4.0, 4.25) 

A = $ABS (B) = (4.0, 4.25) 

(2) B = (-2.0, -1.0) 

A - $ABS (B) - (1.0, 2.0) 


(3) B = (-1.5, 4.0) 

A = $ABS (B) - (0, MAX (-1.5, 4.0)) 
= (0, 4.0) 


12. A = $SCH (B, C) A: = (MAX((CL, BL), MIN(CR, BR)) Intersection 
Ex. B = (4.0, 5.0) $SCH can produce illegal intervals. This can be 
C = (-1 .5 , 2 .0) proved with ICOR . 

A = $SCH (B, C) = (MAX (4.0, -1.5), MIN (5.0, 2.0) 

= (4.0, 2.0) (- illegal interval) 

B - (1.0, 2.0) 

C = (-1.5, 5.0) 

A = $SCH (B, C) = (1.0, 2.0) 


13. A - $VER (B, C) A: = (MIN (CL, BL), MAX(CR, BR)) Union 

B = (1.5, 2.5) 

C = (1.0, 5.0) 

A = $VER (B, C) = (MIN (1.0, 1.5), MAX (5.0, 2.5)) 

= (1.0, 5.0) 

14. IF (IZER (A)) \ “““ Test to see if zero (0) is contained in A 

If IZER(A) < 0 GO TO" 

If IZER(A) = 0 GO TO ^ 

If IZER(A) > 0 GO TO 
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Test to see: 


* : if AR < 0 (if AR <0,0^ A) 

*** : if AL > 0 (if AL > 0 , 0 ^ A) 

** : else 0 e A (zero contained in A) null interval 

15. IF (ICOR(A)) % Interval Test 

IF ICOR(A) < 0 GO TO" 

IF ICOR(A) = 0 GO TO ** 

IF ICOR(A) > 0 GO TO 

Test to see: 

* : if AL > AR (illegal interval) 

~~ : if AL = AR 

: if AL < AR 

16. $ENT — Test on being contained in 
CALL $ENTC A, B, &“1, &~2, &*3 f & '4) 

“ 1: AL > BL and AR < BR (A is contained in B) 

*2: AL = BL and AR = BR (A and B are equal) 

' 3: AL < BL and AR > BR (B is contained in A) 

~ 4: A and B are not fully contained in one another or are not contained in 
each other at all. 


17. $SUM - Build a sum with a scalar product. 

CALL $SUM(A, I, J, K) where I = J, K 

$SUM store A, I, J, K and sets up a loop made to perform addition. 


Along with the call to $SUM a multiplication statement has to follow 


immediately. 

CALL $SUM (A, I, J, K) 
D = $MUL(B, C) 


Ex: 

A = (0.0, 0.0) 

B(l) = (1.0, 2.0) C(l) = (2.0, 3.0) 
B(2) = (3.0, 4.0) C(2) = (4.0, 5.0) 
B(3) = (5.0, 6.0) C(3) = (6.0, 7.0) 


Example: 

CALL $SUM (A, I, 1, 3) 
D = $MUL (B(I), C(I)) 


J = 1, K = 3 

CALL $SUM (A, I, J, K) 
D = $MUL (B(I), C(I)) 


Correspondent 


D = (0.0, 0.0) + (2.0, 6.0) 

- ( 2 . 0 , 6 . 0 ) + ( 12 . 0 , 20 . 0 ) 

= (14.0, 26.0) + (30.0, 42.0) 
= (44.0, 68.0) 


AKKU = A 
DO 1 I = 1, 3 
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1 AKKU = AKKU + $MUL (B(I), C(I)) 
D = AKKU 


D, AKKU, and A are intervals; B, C are interval arrays. The sequence of $SUM 
and $MUL cannot contain interval functions. 

18. $INP A = $INP (’ (INPUT >’) conversion routine 

(INPUT) is an interval number, that is produced from an alphanumeric field 
of input values surrounded by quotation marks. The routine creates an 
interval from the set of numbers in the alphanumeric field by choosing 
the smallest number for the left endpoint and the largest number for the 
right endpoint. The interval number is greater than or equal to 1 (> 1) 
in absolute value. The real or integer numbers are separated by commas 
but can contain no blanks. The beginning and end of the interval number 
must contain at least one blank. The left side of the output is the minimum 
of the real or integer numbers. The right side of the output contains the 
maximum of the real or integer numbers. 

Example: 


A = $INP (’# 1.3, 1.56#') -> 
A = $INP ('# -2., 4#') 

A = $INP (’!#') 


A = (1.3, 1.56) 

A - (-2., 4.) 

No answer a blank doesn’t preceed "1" on 
left side 


A = $INP ('# 1E45, 23, 2.5#') _ A = (2.5, 1.E45) 


The arithmetic functions can produce exponent overflow and divide check 
interrupts. Certain interval functions can produce illegal intervals and only two 
of the subroutines in this package can determine if the output is illegal, IZER 
and ICOR. 


In all the arithmetic interval functions the input is tested to determine if 
it is illegal. No error messages are generated just a return to the calling 
program. 


DEFINITION OF (INPUT): 


(INPUT) 

: = “ (INTER) ~ 

(INTER) 

: = (REAL) 


: = < INTER >, (REAL) 

REAL 

: = (FLOAT) 


: = (FLOAT) (EX) < INTEG) 
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<FLOAT> 

: = <INTEG> 


: = <INTEG> • <ZIFF> 

(INTEG) 

:= (SIGN) (ZIFF) 

< ZIFF) 

: = (DIGIT) 


:= (ZIFF) (DIGIT) 

< DIGIT) 

: = 1234567890 

<SIGN> 

: = + - < > 

<EX> 

: = E 


: = D 

* 

At least one blank 

< > 

The empty string. 


V. MATHEMATICS USED IN PACKAGE 

If * is one of the symbols +, ■>, /, the arithmetic operations of intervals 

are defined by: 

[a, b] ■ * [c, d]={x*y|a<x<b, c < y < d} 

[a, b] + [c, d] = [a + c, b + d] 


[ a, b] ■ - [c, d] = [a - d, b - c] 


[a, b] • [c, d] = [min (ac, ad, be, bd), max (ac, ad, be, bd)] 


and if 0 ^ [c, d] then 


[a, b] ■/ [c, d] = [ a, b] 
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Normalized floating point interval arithmetic is used in this package. Under 
IBM O/S 360 a floating point number is represented as 


bit 

0 

1 - 7 

8 63 


sign 

characteristic 

fraction (binary) 



exponent 

mantissa 


where bit = 0 sign of number (0-positive, 1 negative) 
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1 - 7 = the exponent indicating the power of 16 by which the fraction 
is multiplied. A zero in bit 1 indicates negative exponent, 

1 indicates positive exponent. It is coded in excess 64- 
notation (decimal). That is, subtract decimal equivalent 
of characteristic by 64 to obtain the actual characteristic . 
e.g. 0 1 1000000 1 = 64 
64 - 64 = o = 16 

8 - 63 = fraction or mantissa. The binary points of the fraction is 

just before bit position 8. That is to say that a 1 in bit posi- 
tion 8 represents 2" 1 , a 1 in bit position 9 represents 2“ 2 , etc 

The magnitudes of a floating point number is 16 -65 (5.4 x 10 -79 ) 

(1 - 16" 14 ) • 16 63 (7.2 x 10 7S ) 

Example: 

1 1000001 01010100 00000000 00000000 00000000 etc. 

Sign Characteristic 

(-1) 65 

65 - 64 = 1 

- 1 x 16 1 x (2“ 2 + 2“ 4 + 2~ 6 ) = -2 4 (2 -2 + 2" 4 + 2" 6 ) 

= - (2 2 +2° + 2“ 2 ) 

= - (4 + 1 + 0.25) = -5.25 

Normalized floating point numbers have a non-zero high order hexadecimal 
fraction digit. One or more high order hexadecimal fraction digits which are 
zero makes the number unnormalized. Normalization consists of shifting the 
fraction left until the high order hexadecimal digit is non-zero and reducing 
the characteristic by the number of hexadecimal digits shifted. For multiplica- 
tion and division the operands are normalized before the operations are per- 
formed. For the remaining operands (+, -) the result is normalized. 

A = 436D1000 00000000 = normalized number 

B = 4306D100 00000000 = unnormalized 

= 426D1000 0000000 after normalization. 
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VI. RESTRICTIONS 


1. Only basic arithmetic operations and basic set theory applications are per- 
formed. No trigonometric operations are performed. 

2. Uses long (64 bits) floating point interval arithmetic. 

3. The arithmetic performed can produce illegal intervals. 

4. The arithmetic performed can cause arithmetic overflow and divide check 
interrupts. 
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-APPENBKE 


LEVEL 20.1 (AUG 71 ) 


OS/ 360 FORTRAN H 


DATE 72.108/09. 51.35 


COMPILER OPTIONS - NAME = MAI N , OP T=0 1 , LI NEC NT = 58 , S I ZE = OOOOK , 

SOURCE, EBCDIC, NOLI ST, NODECK, LOAD, MAP, NOEDIT, I D, XREF 

ISN 0002 IMPLIC I T COMPL EX *16 ( A-H , $ ) ,RE A L*8 (0-Z) ,_I NT EGER ( I -N) 

ISN 0003 DIMENSiON C ( 20 ) , D ( 20 ) , E ( 20 ) 

ISN 0004 WRITE ( 6, 101 ) 

ISN 0005 -1=1 

ISN 0006 P = 0 ♦ 5D0 

ISN 0007 A= ( 2 • ODO , 5 .ODO ) '' 

ISN 0008 B= ( 1.000,8. ODO) 


ISN 0009 N=3 

ISN 0010 V = DL I ( A ) 


ISN 

0011 

S=DRE ( A ) 

ISN 

0012 

T=DL I ( R ) 

ISN 

0013 

U=DR E ( 8 ) 

ISN 

0014 

WRITE! 6,203) V,S,T,U,N,P 

ISN 

0015 

WR I TE ( 6, 222) 

ISN 

0016 

D( 1 ) = $ I T$ ( N ) 

ISN 

0017 

WRITEI6, 108) 

ISN 

0018 

R=DL I (D(I)) 

ISM 

0019 

S =DR E ( D ( I ) ) 

ISN 

0020 

WRITE! 6i.2002) I , R, S 

ISN 

0021 

1=1 + 1 

ISN 

0022 

DI2)=$DT$(P) 

ISN 

0023 

WRITEI6, 109) 

ISN 

0024 

R=DL I (0(1) ) 

ISN 

0025 

S=DRE(D( I ) ) 

ISN 

0026 

WRI TE ( 6,2002 ) I,R,S 

ISN 

0027 

1=1 + 1 

ISN 

0028 

D ( 3 ) = $ ADD ( A , B ) 

ISN 

0029 

WRITE) 6, 102) 

ISN 

0030 

R = DL I (0(D) 

ISN 

0031 

S =DRE ( D ( I ) ) 

ISN 

0032 

WRITE ( 6,2002 ) I,R,S 

ISN 

0033 

1=1 + 1 

ISN 

0034 

D ( 4 ) = $ SUB ( A , B ) 

ISN 

0035 

WRITE! 6, 103) 

ISN 

0036 

R=DL I (0(1)) 

ISN 

0037 

S=DRE(D< I ) ) 

ISN 

0038 

WRITE! 6,2002) I,R,S 

ISN 

0039 

1=1 + 1 

ISN 

0040 

D ( 5 ) = $MUL ( A , B ) 

ISN 

0041 

WRITEI6, 104) 

ISN 

0042 

R=DL I ( D ( I ) ) 

ISN 

0043 

S = DRE ( D ( I ) ) 

ISM 

0044 

WR I TE i 6,2002 ) i . R , S 

ISN 

0045 

1=1 + 1 

ISN 

0046 

D( 6) =$D I V ( A , B ) 

ISN 

0047 

WRITE! 6, 105) 

ISN 

0048 

R = DL I ( D ( I ) ) 

ISN 

0049 

S=DRE(D< I ) ) 

ISN 

0050 

WRITE! 6,2002) I.R.S 

ISN 

0051 

WRITE(6,223) 

ISN 

0052 

1 = 1 

ISN 

0053 

C(1)=SP0T(A,N) 

ISN 

0054 

WR I TE ( 6, 1 10 ) 

ISN 

0055 

R =DL I ( C ( I ) ) 
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ISN 

0056 

S=DRE ! C( I ) ) 

ISN 

0057 

WRITE! 6,2002) I , R, S 

ISN 

0058 

1=1 + 1 

ISN 

0059 

C(2)=$P0T!B,N> 

ISN 

00 60 

WR I TE ( 6, 111) 

ISM 

0061 

R=DL I ( C ( I ) ) 

ISN 

0062 

S=DRE ( C ( I ) ) 

I SN 

0063 

WRITE! 6,2002) I,R»S 

ISN 

00 64 

1=1 + 1 

ISN 

0065 

C ( 3 ) = $M IN ( A ) 

ISN 

00 66 

WRI TE ( 6, 1 1 2 ) 

ISN 

00 67 

R=DL I ( C ! I ) ) 

ISN 

00 68 

S=DR E ( C ( I ) ) 

ISN 

00 69 

WRITE ( 6,2002 ) I,R,S 

ISN 

0070 

1=1 + 1 

ISN 

0071 

C(4)=$MIN(B) 

ISN 

0072 

WR I TE ( 6, 113) 

I SN 

0073 

R=DL I ( C ( I ) ) 

ISN 

0074 

S=DRE ( C ( I ) ) 

I SN 

0075 

WRITE ( 6,2002) I,R,S 

ISN 

0076 

1=1 + 1 

ISN 

0077 

C( 5)=$ABS( A) 

ISN 

0078 

WRI TE ( 6, 114) 

ISN 

0079 

R=DL I ( C ( I ) ) 

ISN 

0080 

S =DRE ( C ( I ) ) 

ISN 

0081 

WRITE 1 6,2002) I,R»S 

ISN 

0082 

1=1 + 1 

ISN 

0083 

C< 6) = $ABS ( B ) 

ISN 

0084 

WRITE! 6, 115) 

ISN 

0085 

R=DL I ( C ( I ) ) 

ISN 

0086 

S=DRE!C( I ) ) 

ISM 

0087 

WRITE ( 6,2002) I,R,S 

ISN 

0088 

1=1 + 1 

ISN 

0089 

C ( 7 ) = S SCH ( A , B ) 

ISN 

0090 

WRITE (6, 107) 

ISN 

0091 

R = OL I < C ( I ) ) 

ISN 

0092 

S=DRE(C( I ) ) 

ISM 

0093 

WRITE! 6,2002 ) I,R,S 

ISM 

0094 

I F ( I COR ( C ( 7 ) ) > 10,20,20 

ISN 

0095 10 IF(R.GT.S) WRITE<6,190> R,S 

ISM 

0097 

1=1 + 1 

ISN 

0098 20 C ( 8 ) = $ VER ( A * B ) 

ISN 

0099 

WRITEI6, 106) 

ISN 

0100 

R = DL I < C ( I ) ) 

ISN 

0101 

S=DRE ( C ( I ) ) 

ISN 

0102 

WRITE! 6,2002) I,R,S 

ISM 

0103 

I F ( I ZER ( I ) ) 28,30,29 

ISN 

0104 28 IF (S.LT.O) GO TO 35 

ISN 

0106 29 IF(R.GT.O) GO TO 35 

I SN 

010-8 30 WRITE! 6*210) 


C TEST ON BEING CONTAINED IN 

ISN 

0109 35 WR I TE ( 6, 2 1 7 ) 

ISN 

0110 

CALL SENT! A, B, 636, £37, £38, 639) 

ISN 

0111 36 WRITE! 6,211 )' 

ISN 

0112 

GO TO 40 

ISN 

0113 37 WR I TE ( 6, 2 1 2 ) 
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ISM 0114 
ISM 0115 
ISM 0116 
ISM 0117 
ISM 0118 
ISM 0119 
ISM 01 20 
ISM 0121 


00 TO 40 

38 WRITE <6, 21 3) 

00 TO 40 

39 WRITE (6*214) 

40_ R= ( 0.000,0.000) 

WRITE (6, 101) 

T=nti(B) 

U=DRE < B> 


ISM 0122 
ISM 0123 
ISM 0124 
ISN 0125 
ISM 0126 
ISM 0127 
ISN 0128 
ISN 0129 
ISN 0130 
ISN 0131 
ISM 0132 
ISN 0133 
ISN 0134 
ISN 0135 
ISN 0136 


WR I TE ( 6, 203 ) V,S,T,U,N,p 
WRITE ( 6, 217 ) 

CALL SENT ( A , B » E41, £42 * f.43* £44 ) 
41 WRITE(6,211) 

GO Til 44 

'42'" WR I TF ( 6, 21 2 ) 

GO TO 45 

43 WRITE(6*213) 

GO TO 45 

44 WRITE ( 6, 214) 

45 CONTINUE 

CALL SSUM) B, I , 1, 5 ) 

F = SMUL (C(I) ,D(I ) ) 

R = 0L I (F) 

S = DRE ( F ) 


ISN 0137 WRITE(6*218) 

ISN 0138 WR I TE ( 6* 2 1 5 ) R*S 

C CONVERSION ROUTINE 

ISN 0139 _ C (..!_) = SI.NP ( 1 1.3,1.56 •) 

ISN 0140 C ( 2 ) =$ I MP ( • -2. ,4 1 ) 

ISN 0141 C ( 3 ) = S I NP ( * 1 ') 

ISN 0142 C ( 4 ) = $ I M P ( • 1 E45 , 2 3,2.5 ') 

ISM 0143 CJ_5i=$DTS(5.5D0)_ 

ISN 0144 C ( 6 ) = $ I TS ( 9 ) 


ISN 

0145 


WRITE (6, 216) 

ISM 

0146 


DO 9 1=1,6 

I SN 

0147 


R=OL I ( C ( I ) ) 

ISM 

0148 


S =DR E ( C ( I ) ) 

I SN 

0149 

9 

WRITE ( 6,2002) I,R,S 

ISN 

0150 

2002 

FORMAT) IX , I 5 , • LEFT= ',D24.16,' R I GHT = ',024.16) 

ISM 

0151 

10 1 

FORMAT! ' 1INTERVAL ARITHMETIC PACKAGE') 

I SN 

0152 

102 

FORMAT) ' OA+B' ) 

I SN 

0153 

10 3 

FORMAT) 'OA-R' ) 

I SN 

0154 

104 

FORMAT) ' OA*B ' ) 

I SN 

0155 

10 5 

FORMAT ( 'OA/B' ) 

ISM 

0156 

106 

FORMAT) ' OA UNION B ' ) 

ISN 

015.7 

10 7 

FORMAT) 'OA INTERSECT B') 

TQM 

01 *5 fl 

i no 

cnpMATMntmfMt n 

ISM 

0159 

109 

FORMAT) 1 OSDTS ( f>) ' ) 

ISN 

0160 

1 10 

FORMAT) ' OA**N' ) 

ISN 

0161 

111 

FORMAT) ' OB**N ' ) 

I SN 

0162 

112 

FORMAT) ' OSMIN) A) ' ) 

ISN 

0163 

113 

FORMAT) 'OSMIN(B) ' ) 

ISN 

0164 

114 

FORMAT) ' OSABS) A) ' ) 

ISN 

0165 

115 

FORMAT) 'OSABS (B) ' ) 

ISN 

01 66 

190 

FORMAT)' THIS IS AN ILLEGAL INTERVAL ( ' »G 24 • 1 6 , ', ',024.16, ' ) ' 

I SN 0167 

...20.3 . 

.. FORMAT)!/, ' A= ' ,024.16, ' t ' ,D24.16/1X, 'B= ’ ,024. 16 , ' , • ,024. 16/ 




1 IX, 'N= ' , I 5 / 1 X * ' P = ',024.16) 
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ISN 

0168 

210 

FORMAT! 

ZERO IS CONTAINED IN THE FOLLOWING 

INTERVAL ( ' ,G24.16, ' , ' 




1G24.16, 

TH 


ISN 

0169 

211 

FORMAT! 

A IS CONTAINED IN B' ) 


ISM 

0170 


FORMAT! 

A AND B ARE EQUAL ' ) 


ISN 

0171 

ipWCT-' 

FORMAT! 

B IS CONTAINED IN A' ) 


ISN 

0172 


FORMAT! 

A AND B ARE NOT FULLY CONTAINED IN 

EACH OTHER ' ) 

ISN 

0173 


FORMAT! 

LEFT= • , D24. 16, 2X, ' RIGHT= ',024.16) 


ISN 

0174 

■Ha 

FORMAT! 

0$ I N P CONVERSION ROUTINE') 


ISN 

0175 

mu 

FORMAT! 

0$EN T-TES T ON BEING CONTAINED IN' ) 


ISN 

0176 

218 

FORMAT! 

OSUM PRODUCED WITH SCALAR PRODUct') 


ISN 

0177 

222 

FORMAT! 

OVALUES FOR ARRAY D( I) ' ) 


ISN 

0178 

22 3 

FORMAT ( 

OVALUES FOR ARRAY C (I ) ' ) 


ISN 

0179 


STOP 



ISN 

0180 


END 
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INTERVAL ARITHMETIC PACKAGE 


A= 0 • 2 0000 000 00 000 0000 01, 0.8000000000000000D 01 

B= 0.0 ♦ 0.0 

N= 3 

P= 0.5000000 00 00 0000 OD 00 


SENT-TEST ON BEING CONTAINED IN 


A AND B ARE NOT FULLY CONTAINED IN EACH 

OTHER 

SUM PRODUCED WITH SCALAR PRODUCT 


LEFT= -0. 6850000000000000D 02 RIGHT= 

0.8730000000000000D 03 

$ I NP CONVERSION ROUTINE 



1LEFT= 0.1300000000000000D 01 RIGHT= 0 . 1 560000000006500D 01 
2LEFT= -0. 200 0000000 00 0000 D 01 RIGHT= 0.40000 0000000000QD 0 1 
3LEFT= -0.2000000000000000D 01 RIGHT= 0 . 400 00000 0000 OO'OOb 6T 
4LEFT= 0. 2 50 00 00000 00 0000 D 01 R I GHT= 0 ♦ 100 0000000000 000 D 46 
5L E FT= 0 . 55 OOOOOOOOOOOOOOD 01 R I GHT= 0.5500000b0000b600’6 01 
6LEFT= 0.90000000000000000 01 R I GH T= 0 .900 00000000000000 01 
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INTERVAL ARITHMET_IC_ PACKAGE 

A= 0.20000000000000000 01, 0 . 5000000000000000D 01 

_B= 0 . lOOOOOOOOO OOOOOOD 01, 0.8 000000 OOOOO OOOO D 0 1 1 


N = 3 

P_= 0 .500000000000 00000 00 


VALUES FOR 

ARRAY D ( I ) 





$IT$(i\l) 

1LEFT= 

0.30000000000000000 

01 

R I GHT= 

0. 30 00 000 OOOOOOOOOD 

01 

$DT$ ( R ) 
2_LEFT= 

0.50000000000000000 

.00.. 

R I GHT= 

0.5000000000000000D 

op 

A+B 






3L EFT = 

0 .3000000000000000D 

01 

R I GHT= 

0.13000000000000000 

02 

A-R 

4LEFT.= . 

-0.60000000000000000 

.01 

R I GHT = 

.. ^0. 40000000.00000000 D 

.0 L. 

A*B 






5LEFT= 

0.20000000000000000 

01 

R I GHT = 

0.40000000000000000 

02 

A/R 

6L E FT = 

0.2 5000000000000000 

„00 

R I GHT= 

0.5000000000000000 0 

01 

VALUES FOR 

ARRAY C(I) 





A * * N 






1 L EFT = 

0 .800 00000000000000 

01 

R I GHT = 

0. 12 50 00000000000 OD 

03 

R**N 

21. F FT = 

o. i oooooonqoooooooo 

01 

R I G H T = 

0.51200000000000000 

03 

$MIM( A ) 

3L EFT = 

-0 . 500 00 000 00 0 0 0 0 00 D 

01 

R I GHT = 

-0.20000000000000000 

01 

SMIN( R) 

4LEFT= 

- 0 .80000 00 poo 00 0.0 0 0 0 

01 

R I.G H T= 

- Q.100 00.0 00000 00000 0 

.0.1 

SABS (A ) 






5L EFT = 

0.20000000000000000 

01 

R I GH T = 

0.50000000000000000 

01 

$ A B S ( B ) 

6L E FT = 

0.10000000000000000 

01 

R I GHT= 

0.80000000000000000 

01 

A INTERSECT 

R 





7LEFT= 

0.20000000000000000 

01 

R I GHT = 

0.5000000000000000D 

01 

A UNION B 
8L E FT = 

0. lOOOOOOOOOOOOOOOD 

01 

R I GHT= 

0.80000000000000000 

01 

$_E NT-TEST ON BEING CONTAI NED .. IN . .. 

• 





A IS contained in r 
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